A Transit Timing Analysis of Mine RISE Light Curves of the Exoplanet System TrES-3

We present nine newly observed transits of TrES-3, taken as part of a transit timing program using the RISE instrument on the Liverpool Telescope. A Markov-Chain Monte Carlo analysis was used to determine the planet star radius ratio and inclination of the system, which were found to be R-p/R-star = 0.1664(-0.0018)(+0.0011) and i = 81.73(-0.04)(+0.13), respectively, consistent with previous results. The central transit times and uncertainties were also calculated, using a residual-permutation algorithm as an independent check on the errors. A re-analysis of eight previously published TrES-3 light curves was conducted to determine the transit times and uncertainties using consistent techniques. Whilst the transit times were not found to be in agreement with a linear ephemeris, giving chi(2) = 35.07 for 15 degrees of freedom, we interpret this to be the result of systematics in the light curves rather than a real transit timing variation. This is because the light curves that show the largest deviation from a constant period either have relatively little out-of-transit coverage or have clear systematics. A new ephemeris was calculated using the transit times and was found to be T-c(0) = 2454632.62610 +/- 0.00006 HJD and P = 1.3061864 +/- 0.0000005 days. The transit times were then used to place upper mass limits as a function of the period ratio of a potential perturbing planet, showing that our data are sufficiently sensitive to have probed sub-Earth mass planets in both interior and exterior 2:1 resonances, assuming that the additional planet is in an initially circular orbit.

A transit timing analysis of seven RISE light curves of the exoplanet system HAT-P-3

We present seven light curves of the exoplanet system HAT-P-3, taken as part of a transit timing programme using the rapid imager to search for exoplanets instrument on the Liverpool Telescope. The light curves are analysed using a Markov chain Monte Carlo algorithm to update the parameters of the system. The inclination is found to be i = 86.75+0.22-0.21°, the planet-star radius ratio to be Rp/R* = 0.1098+0.0010-0.0012 and the stellar radius to be R* = 0.834+0.018-0.026Rsolar, consistent with previous results but with a significant improvement in the precision. Central transit times and uncertainties for each light curve are also determined, and a residual permutation algorithm is used as an independent check on the errors. The transit times are found to be consistent with a linear ephemeris, and a new ephemeris is calculated as Tc(0) = 2454856.70118 +/- 0.00018 HJD and P = 2.899738 +/- 0.000007 d. Model timing residuals are fitted to the measured timing residuals to place upper mass limits for a hypothetical perturbing planet as a function of the period ratio. These show that we have probed for planets with masses as low as 0.33 and 1.81 M? in the interior and exterior 2:1 resonances, respectively, assuming the planets are initially in circular orbits.

Modelling the mechanical response of the Reed-mouthpiece-lip system of a clarinet. Part I. A one-dimensional distributed model

The motion of a clarinet reed that is clamped to a mouthpiece and supported by a lip is simulated in the time-domain using finite difference methods. The reed is modelled as a bar with non-uniform cross section, and is described using a one-dimensional, fourth-order partial differential equation. The interactions with the mouthpiece Jay and the player's lip are taken into account by incorporating conditional contact forces in the bar equation. The model is completed by clamped-free boundary conditions for the reed. An implicit finite difference method is used for discretising the system, and values for the physical parameters are chosen both from laboratory measurements and by accurate tuning of the numerical simulations. The accuracy of the numerical system is assessed through analysis of frequency warping effects and of resonance estimation. Finally, the mechanical properties of the system are studied by analysing its response to external driving forces. In particular, the effects of reed curling are investigated.

Plug Force Monitoring for the Control and Optimisation of the Thermoforming Process

Thermoforming processes generally employ sheet temperature monitoring as the primary means of process control. In this paper the development of an alternative system that monitors plug force is described. Tests using a prototype device have shown that the force record over a forming cycle creates a unique map of the process operation. Key process features such as the sheet modulus, sheet sag and the timing of the process stages may be readily observed, and the effects of changes in all of the major processing parameters are easily distinguished. Continuous, cycle-to-cycle tests show that the output is consistent and repeatable over a longer time frame, providing the opportunity for development of an on-line process control system. Further testing of the system is proposed.

Escape times for rigid Brownian rotators in a bistable potential from the time evolution of the Green function and the characteristic time of the probability evolution

The greatest relaxation time for an assembly of three- dimensional rigid rotators in an axially symmetric bistable potential is obtained exactly in terms of continued fractions as a sum of the zero frequency decay functions (averages of the Legendre polynomials) of the system. This is accomplished by studying the entire time evolution of the Green function (transition probability) by expanding the time dependent distribution as a Fourier series and proceeding to the zero frequency limit of the Laplace transform of that distribution. The procedure is entirely analogous to the calculation of the characteristic time of the probability evolution (the integral of the configuration space probability density function with respect to the position co-ordinate) for a particle undergoing translational diffusion in a potential; a concept originally used by Malakhov and Pankratov (Physica A 229 (1996) 109). This procedure allowed them to obtain exact solutions of the Kramers one-dimensional translational escape rate problem for piecewise parabolic potentials. The solution was accomplished by posing the problem in terms of the appropriate Sturm-Liouville equation which could be solved in terms of the parabolic cylinder functions. The method (as applied to rotational problems and posed in terms of recurrence relations for the decay functions, i.e., the Brinkman approach c.f. Blomberg, Physica A 86 (1977) 49, as opposed to the Sturm-Liouville one) demonstrates clearly that the greatest relaxation time unlike the integral relaxation time which is governed by a single decay function (albeit coupled to all the others in non-linear fashion via the underlying recurrence relation) is governed by a sum of decay functions. The method is easily generalized to multidimensional state spaces by matrix continued fraction methods allowing one to treat non-axially symmetric potentials, where the distribution function is governed by two state variables. (C) 2001 Elsevier Science B.V. All rights reserved.

Hydrochemical and isotopic tracing of mixing dynamics and water quality evolution under pumping conditions in the mine shaft of the abandoned Frances Colliery, Scotland.

Since 1995, when pumps were withdrawn from deep mines in East Fife (Scotland), mine waters have been rebounding throughout the coalfield. Recently, it has become necessary to pump and treat these waters to prevent their uncontrolled emergence at the surface. However, even relatively shallow pumping to surface treatment lagoons of the initially chemically-stratified mine water from a shaft in the coastal Frances Colliery during two dynamic step-drawdown tests to establish the hydraulic characteristics of the system resulted in rapid breakdown of the stratification within 24 h and a poor pumped water quality with high dissolved Fe loading. Further, data are presented here of hydrochemical and isotopic sampling of the extended pump testing lasting up to several weeks. The use in particular of the environmental isotopes d18O, d2H, d34S, 3H, 13C and 14C alongside hydrochemical and hydraulic pump test data allowed characterisation of the Frances system dynamics, mixing patterns and water quality sources feeding into this mineshaft under continuously pumped conditions. The pumped water quality reflects three significant components of mixing: shallow freshwater, seawater, and leakage from the surface treatment lagoons. In spite of the early impact of recirculating lagoon waters on the hydrochemistries, the highest Fe loadings in the longer-term pumped waters are identified with a mixed freshwater–seawater component affected by pyrite oxidation/melanterite dissolution in the subsurface system.

ULTRACAM photometry of the eclipsing cataclysmic variable OU Vir

We present high-speed, three-colour photometry of the faint eclipsing cataclysmic variable OU Vir. For the first time in OU Vir, separate eclipses of the white dwarf and the bright spot have been observed. We use timings of these eclipses to derive a purely photometric model of the system, obtaining a mass ratio of q=0.175+/-0.025, an inclination of i=79.degrees2+/-0.degrees7 and a disc radius of R-d/a=0.2315+/-0.0150. We separate the white dwarf eclipse from the light curve and, by fitting a blackbody spectrum to its flux in each passband, obtain a white dwarf temperature of T=13900+/-600 K and a distance of D=51+/-17 pc. Assuming that the primary obeys the Nauenberg mass-radius relation for white dwarfs and allowing for temperature effects, we also find a primary mass M-w/M-circle dot=0.89+/-0.20, a primary radius R-w/R-circle dot=0.0097+/-0.0031 and an orbital separation a/R-circle dot=0.74+/-0.05.

The Exotic Eclipsing Nucleus of the Ring Planetary Nebula SuWt 2

SuWt 2 is a planetary nebula (PN) consisting of a bright ionized thin ring seen nearly edge-on, with much fainter bipolar lobes extending perpendicularly to the ring. It has a bright (12th magnitude) central star, too cool to ionize the PN, which we discovered in the early 1990s to be an eclipsing binary. Although it was anticipated that there would also be an optically faint, hot, ionizing star in the system, a spectrum from the International Ultraviolet Explorer (IUE) did not reveal a UV source. We present extensive ground-based photometry and spectroscopy of the central binary collected over the ensuing two decades, resulting in the determination that the orbital period of the eclipsing pair is 4.9 days, and that it consists of two nearly identical A1 V stars, each of mass ~2.7 M sun. The physical parameters of the A stars, combined with evolutionary tracks, show that both are in the short-lived "blue-hook" evolutionary phase that occurs between the main sequence and the Hertzsprung gap, and that the age of the system is about 520 Myr. One puzzle is that the stars' rotational velocities are different from each other, and considerably slower than synchronous with the orbital period. It is possible that the center-of-mass velocity of the eclipsing pair is varying with time, suggesting that there is an unseen third orbiting body in the system. We propose a scenario in which the system began as a hierarchical triple, consisting of a ~2.9 M sun star orbiting the close pair of A stars. Upon reaching the asymptotic giant branch stage, the primary engulfed the pair into a common envelope, leading to a rapid contraction of the orbit and catastrophic ejection of the envelope into the orbital plane. In this picture, the exposed core of the initial primary is now a white dwarf of ~0.7 M sun, orbiting the eclipsing pair, which has already cooled below the detectability possible by IUE at our derived distance of 2.3 kpc and a reddening of E(B - V) = 0.40. The SuWt 2 system may be destined to perish as a Type Ia supernova.

EXPERIMENTAL STUDIES OF THE MNF PHOTORESONANT EXCITATION SOFT-XRL SCHEME

Photoresonant excitation has been examined as a mechanism for enhancing the maximum gain achievable on a H-like recombination scheme. Modelling of the system, using the MED101/NIMP package, has shown that the photoresonant excitation of the 1s-3p (Ly(beta)) transition at 12.644 angstrom in a H-like fluorine plasma by a Be-like manganese transition at 12.643 angstrom offers the possibility of considerable gain enhancement on the 3d-2p Balmer alpha at 80.91 angstrom over the simple H-like system. Some experimental evidence for a transitory increase in intensity on the Balmer-alpha has been obtained for a fluorine-coated C-fiber system pumped by irradiating a manganese stripe in close proximity to the fiber.

Dynamics of the entanglement spectrum in spin chains

We study the dynamics of the entanglement spectrum, that is the time evolution of the eigenvalues of the reduced density matrices after a bipartition of a one-dimensional spin chain. Starting from the ground state of an initial Hamiltonian, the state of the system is evolved in time with a new Hamiltonian. We consider both instantaneous and quasi adiabatic quenches of the system Hamiltonian across a quantum phase transition. We analyse the Ising model that can be exactly solved and the XXZ for which we employ the time-dependent density matrix renormalisation group algorithm. Our results show once more a connection between the Schmidt gap, i.e. the difference of the two largest eigenvalues of the reduced density matrix and order parameters, in this case the spontaneous magnetisation.

Full characterization of the quantum linear-zigzag transition in atomic chains

A string of repulsively interacting particles exhibits a phase transition to a zigzag structure, by reducing the transverse trap potential or the interparticle distance. Based on the emergent symmetry Z2 it has been argued that this instability is a quantum phase transition, which can be mapped to an Ising model in transverse field. An extensive Density Matrix Renormalization Group analysis is performed, resulting in an high-precision evaluation of the critical exponents and of the central charge of the system, confirming that the quantum linear-zigzag transition belongs to the critical Ising model universality class. Quantum corrections to the classical phase diagram are computed, and the range of experimental parameters where quantum effects play a role is provided. These results show that structural instabilities of one-dimensional interacting atomic arrays can simulate quantum critical phenomena typical of ferromagnetic systems.

Optimizing the efficiency of the sub-map technique for large-scale simultaneous localization and mapping

A technique for optimizing the efficiency of the sub-map method for large-scale simultaneous localization and mapping (SLAM) is proposed. It optimizes the benefits of the sub-map technique to improve the accuracy and consistency of an extended Kalman filter (EKF)-based SLAM. Error models were developed and engaged to investigate some of the outstanding issues in employing the sub-map technique in SLAM. Such issues include the size (distance) of an optimal sub-map, the acceptable error effect caused by the process noise covariance on the predictions and estimations made within a sub-map, when to terminate an existing sub-map and start a new one and the magnitude of the process noise covariance that could produce such an effect. Numerical results obtained from the study and an error-correcting process were engaged to optimize the accuracy and convergence of the Invariant Information Local Sub-map Filter previously proposed. Applying this technique to the EKF-based SLAM algorithm (a) reduces the computational burden of maintaining the global map estimates and (b) simplifies transformation complexities and data association ambiguities usually experienced in fusing sub-maps together. A Monte Carlo analysis of the system is presented as a means of demonstrating the consistency and efficacy of the proposed technique.

The optimized geometric structure of the (0001) surface of α-Fe 2O3

A tridimensional model of α-Fe₂O₃ and models of (0001) and (1102) surfaces on it were built. Then the structural optimization of the (0001) surface was presented which explored the influence of the system scale and the terminal surface configuration. Four different models including two different system scale structures (MODEL□ and MODEL□) and two different terminal structures (MODEL□ and MODEL□) were analyzed in this paper. It was concluded that the boundary effect was more important in a smaller system in the structure optimization. And the Fe-terminated was more stable than the O-terminated structure which was agreed with the experiences, this structural model can be used in further work including the monatomic adsorption/desorption and the chemical reactions on this surface.